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Introduction 

Phase  contrast  imaging  differs  from  conventional  x-ray  shadowgraphy  in  the 
mechanism  of  contrast  generation:  while  conventional  shadowgraphy  depends  on 
absorption  of  x-rays,  phase  contrast  imaging  is  based  on  phase  changes  as  x-rays  traverse 
a  body  resulting  in  wave  interference  that  result  in  intensity  changes  in  the  image.  The 
complex  index  of  refraction  n  for  any  material  is  given  by 

n  =  1— 5— i(3  (1) 

where  the  imaginary  part  /?of  the  index  describes  x-ray  absorption  and  the  real  part  8 
describes  the  phase  shift  undergone  by  the  x-radiation  as  it  passes  through  a  material. 
Fresnel-Huygens  theory  '"3,  which  governs  image  formation  in  general,  describes  contrast 
in  an  image  as  dependent  on  both  absorption  and  phase  contrast.  However,  only  when 
sources  with  high  spatial  coherence,  such  as  synchrotrons  4'6,  microfocus  x-ray  tubes 7,  or 
laser  plasma  x-ray  sources  8,9are  employed  is  the  phase  contrast  component  of  the  image 
visible.  Coherent  sources  give  an  image  of  a  dense  object  embedded  in  lower  density,  soft 
tissue  as  a  shadow  of  the  object,  caused  by  absorption,  surrounded  at  its  perimeter  by 
light  and  dark  interference  fringes  arising  from  rapid  phase  variations  in  the  radiation  at 
the  interface  between  the  two  media. 

A  detailed  analysis  of  the  relative  contrast  produced  in  an  image  for  a  fixed 
difference  in  density  shows  that  phase  contrast  is  more  sensitive  than  absorption 
throughout  most  of  the  spectral  region  commonly  used  for  diagnostic  x-ray  imaging10. 

For  instance,  the  thickness  of  water  needed  to  produce  a  1%  phase-contrast  at  36-keV 
photon  energy  is  2500  times  smaller  than  that  required  to  produce  the  same  contrast 
though  absorption10. 

Phase  contrast  in  an  image  can  be  markedly  enhanced  by  recording  two  images  of 
an  object,  one  where  the  object  is  displaced  slightly  through  acoustic  radiation  pressure11, 
and  a  second  where  the  object  is  unmoved,  followed  by  subtraction  of  the  two  images. 
Acoustic  radiation  forces  can  arise  from  either  reflection  of  sound  by  an  object  or  as  a 
result  of  acoustic  impedance  changes  (i.e.  a  variation  in  either  density  or  sound  speed)  or 
by  dissipative  processes  within  the  object  resulting  in  absorption  of  the  ultrasound  energy 
1  .  A  number  of  studies13'18  have  shown  that  soft  tissue  can  be  manipulated  by  ultrasonic 
pressure.  Typically,  a  mechanical  actuator  or  a  beam  of  ultrasound  moves  an  object 
whose  position  is  monitored  in  time  with  conventional  pulse-echo  ultrasonic  imaging.  It 
has  been  shown  that  variations  in  Young’s  modulus,  which,  in  part,  determines  sound 
speed,  permits  acoustic  differentiation  of  tissue.  It  is  known,  for  instance,  that  breast 
tumors  have  a  significantly  different  Young’s  modulus  from  surrounding  tissue,  and  that 
tumors  move  as  rigid  bodies  in  response  to  acoustic  radiation  pressure  14’16'19.  Intense 
ultrasound  is  used  to  exert  forces  on  objects  within  a  body  giving  displacements  on  the 
order  of  tens  to  hundreds  of  microns.  Subtraction  of  images  made  with  and  without  the 
ultrasound  field  gives  an  image  that  removes  low  spatial  frequency  features  and 
highlights  high  frequency  features.  The  method  acts  as  an  acoustic  “contrast  agent”  for 
phase  contrast  x-ray  imaging,  which  in  soft  tissue,  acts  to  highlight  small  density 
changes.  We  have  named  this  method  acoustically  modulated  X-ray  phase  contrast. 
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We  report  here  some  preliminary  experimental  and  theoretical  studies  of  acoustic 
radiation  pressure.  We  have  also  performed  some  additional  proof  of  concept  type 
experiments  combining  ultrasound  and  X-ray  phase  contrast.  Finally,  we  are  also 
developing  the  software  and  hardware  necessary  to  perform  tomographic  studies;  some 
first  theoretical  and  experimental  tests  have  been  completed. 


Body 

Acoustic  Radiation  Pressure 


For  the  purposes  of  the  present  study,  where  reflection  of  the  acoustic  wave  is  considered, 
it  can  be  said  that  for  a  fixed  acoustic  impedance  change,  the  force  exerted  by  an 
ultrasonic  wave  is  proportional  to  the  time  average  of  the  energy  density  of  the  wave  and 
the  area  of  the  object  presented  to  the  sound  field.  An  expression  for  the  acoustic 
radiation  force  on  an  arbitrary  object  has  been  given  by  Westervelt 1 1  and  evaluated  for 
spheres  small  compared  with  the  wavelength  of  the  radiation.  For  spheres  of  any 
diameter,  King  20  has  given  an  expression  that  agrees  with  the  limiting  expression  given 
by  Westervelt.  According  to  King,  the  average  acoustic  radiation  force  exerted  on  a 
sphere  of  radius  a  is  given  by: 


Fac_(ka)2kVp'A(ka) 


(2) 


r 


where  p  and  c  are  the  density  and  speed  of  sound  of  the  medium,  respectively,  p  is  the 
acoustic  pressure,  k  is  the  wave  number  of  the  radiation,  and  A(ka)  is  given  by 
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where  the  functions  Hn  for  small  ka  are  given  by 


n  =  0 : 
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n  =  1 : 

H,  =  4 

f1  +  T“T 

'  (ka)6 

^  2pJ 

n  =  2  : 

H  -  81  ( 

'i+M_ 

2"(kafi 

1+  10 

n  >  2 : 

(ka)' 


= 


(n  +  l)2(fl(2n-l) 

_ Ui _ . 

(ka)4n+2 


(1+(n+lV2n-l)(ka)2+-) 


where  p  is  the  density  of  the  fluid  and  ps  is  the  density  of  the  sphere.  For  ka  >  2,  the 
functions  Hn  are  given  by 


n  =  0: 

H  J  +  W) 


n  =  1  : 


H, 


2(ka)3 


n  >  1: 


where  Jn  is  a  Bessel  function  of  order  n,  and  were  the  argument  of  each  of  the  Bessel 
functions  is  ka. 

Figure  1  is  a  plot  of  force  versus  ka  from  Eq.  2,  with  the  inset  showing  a  portion 
of  the  curve  for  small  ka.  The  parameters  used  in  Eq.  2  were  for  a  Teflon  bead  with  a  = 
1.5  mm,  p  /ps  =  2  ;  and  for  water  with  c  =  1500  m/s  and  p  =  103  kg/m3.  The  curves  in 
Figure  1  show  a  steady  increase  in  the  force  with  increasing  ka 
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(ka) 


Figure  1:  Radiation  force  (in  AO  versus  ka  (dimensionless)  from  Eq.  6.  Inset:  Magnified  view  of  the  first 
part  of  the  curve.  The  portions  of  the  curves  for  ka<2  were  calculated  using  the  small  ka  expressions  for  the 
functions  Hn. 

Experiments  were  carried  out  to  determine  the  deflection  of  a  Teflon  bead  in 
water  using  an  x-ray  imaging  apparatus  to  determine  the  deflection  of  the  bead  as  a 
function  of  acoustic  pressure.  The  x-rays,  generated  by  a  microfocus  x-ray  tube  (Oxford 
Ultrabright  Microfocus  UB-M1),  were  directed  onto  a  sample  cell  consisting  of  a  3-cm 
diameter  PVC  tube  which  had  portions  machined  out  and  replaced  with  Mylar  foil  to 
reduce  x-ray  intensity  losses.  The  x-ray  beam,  as  shown  in  Fig.  2,  propagated  through  a 
He  filled  tube  to  a  Gd202S(Tb)  fiber  optic  scintillation  plate  (Hamamatsu,  Inc.  Model 
J6676)  the  fluorescence  from  which  was  imaged  onto  a  liquid  nitrogen  cooled  CCD 
camera  (Roper  Scientific  Model  7382-0001).  The  intensity  maximum  of  the  emission  x- 
ray  spectrum  was  40  keV.  The  source-to-detector  distance  was  2.6  m,  while  the  sample- 
to-source  distance  was  0.35  m,  giving  a  magnification  of  7.5  and  a  maximum  phase 
contrast  for  spatial  frequencies  corresponding  to  4.3  pm. 
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Figure  2:  Diagram  of  the  experimental  apparatus.  X-radiation  generated  by  a  microfocus  tube  penetrates  a 
sample  and  is  detected  by  a  CCD  camera  that  views  a  scintillation  plate.  The  CCD  camera  is  read  by  the 
computer  which  stores  the  images  and  performs  the  subtractions. 


The  acoustic  transducer  used  to  produce  radiation  pressure  was  a  1  MHz, 
LiNb03-transdueer  located  3  cm  from  the  bead.  One  end  of  the  PVC  sample  tube  was 
terminated  by  a  synthetic,  sound-absorbing  plastic.  The  pulse  train  from  the  function 
generator  (Agilent  Model  33250A)  used  to  drive  the  transducer  consisted  of  15-p.s  bursts 
at  a  repetition  rate  of  133  Hz.  The  signal  from  the  function  generator  was  amplified  by  a 
power  amplifier  with  a  peak  output  power  of  1 .5  kW  and  delivered  to  the  transducer. 
Although  the  low  duty  cycle  of  the  pulse  sequence,  approximately  0.2%,  resulted  in  a 
reduced  time  average  acoustic  power  from  the  transducer,  the  resulting  radiation  pressure 
was  sufficient  to  cause  displacements  that  could  be  easily  imaged  with  the  CCD  camera. 

A  3  mm  diameter  Teflon  bead  was  attached  to  a  glass  fiber  and  suspended  in  the  center  of 
the  sample  cell,  which  was  filled  with  water.  The  glass  fiber  was  calibrated  by  attaching 
it  to  a  precision  translation  stage  and  using  a  laboratory  balance  to  determine  force  versus 
fiber  displacement.  The  data  shown  in  Fig.  3(A)  gave  a  force  constant  of  2.18  x  10~2  N/m 
for  the  fiber.  The  total  force  on  a  bead  suspended  from  the  fiber  is  the  sum  of  the 
restoring  force  from  the  fiber  and  the  force  from  gravity,  which  can  be  written  for  small 
angles  of  deviation  from  the  vertical  as 


p  -  F  + 

exp  glassfiber 


mbead’g-d 


(4) 


where  rribead  is  the  mass  of  the  bead,  g  is  the  gravitational  constant,  L  is  the  length  of  the 
fiber  and  d  is  its  displacement.  For  the  experiments  described  here,  the  mass  of  the  bead 
was  35  mg,  and  the  length  of  the  glass  fiber  was  10  cm.  The  displacement  of  the  bead 
versus  the  time  averaged  electrical  power  delivered  to  the  transducer  is  shown  in  Fig.  3 
(B).  The  break  in  the  slope  of  the  curve  can  be  attributed  to  a  departure  of  the  efficiency 
of  the  transducer  from  linearity  at  high  driving  voltages.  The  acoustic  pressure  generated 


8 
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by  the  transducer  at  various  power  levels  from  the  power  amplifier  was  determined  using 
a  calibrated  LiNiCb  transducer  (LaserSonics  Technologies,  Model  WAT-04).  The  plot  of 
acoustic  pressure  from  the  calibrated  transducer  versus  average  driving  voltage  applied  to 
the  transducer  is  shown  in  Fig.  4(A)  which  gives  a  curve  with  the  same  qualitative 
features  as  the  curve  in  Fig.  3(B),  and  permits  calibration  of  the  transducer  pressure  in 
terms  of  the  power  amplifier  output  driving  voltage.  With  this  calibration,  the  plot  in  Fig. 
3(B)  can  be  converted  to  a  plot  of  force  versus  the  square  of  the  acoustic  pressure  as 
shown  in  Fig.  4(B).  The  theoretical  curve  from  Eq.  6,  also  plotted  in  Fig  4(B),  shows 
reasonably  good  agreement  with  the  experimental  measurements.  The  maximum  force 
recorded  in  the  experiments  was  5  pN. 


-i - - - 1 - « - 1 - » - 1 - > - r- 
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Figure  3:  (A)  Force  on  the  glass  fiber  in  N  versus  displacement  in  m.  (B)  Displacement  of  the  Teflon  bead 
attached  to  the  fiber  versus  average  electrical  power  delivered  to  the  transducer. 
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Figure  4:  (A)  Acoustic  pressure  in  Pa  from  the  calibrated  transducer  versus  driving  voltage  from  the  power 
amplifier.  (B)  Force  on  the  bead  in  N  versus  the  square  of  the  pressure  from  the  transducer. 

Acoustically  Modulated  X-ray  Phase  Contrast  Imaging 

The  method  of  modifying  a  phase  contrast  image  using  acoustic  radiation  pressure 
21,22  consists  of  two  steps:  first,  an  x-ray  image  is  made  with  a  sound  beam  directed  into  a 
body  to  displace  an  object  through  acoustic  radiation  force  and  the  image  stored  in  a 
computer,  second,  another  x-ray  image  is  of  the  object  is  taken,  this  time  without  the 
presence  of  the  sound  beam  and  the  image  is  recorded  in  the  computer.  The  two  images 
are  then  subtracted  pixel  by  pixel  to  give  a  subtracted  phase  contrast  image,  the 
component  of  the  image  from  absorption  contrast  being  largely  eliminated,  leaving  a 
nearly  pure  phase  contrast  image,  inherently  background  and  flatfield  corrected. 

The  principle  of  the  method  was  demonstrated  by  taking  images  of  two  3-mm 
diameter  Teflon  beads  cast  in  a  6.5-cm  long  block  of  agarose.  As  shown  in  Fig.  5  (A),  the 
two  beads  were  both  within  the  field  of  the  x-ray  beam,  but  only  one  bead  was  irradiated 
with  ultrasound.  Figure  5  (B)  shows  an  image  of  the  two  beads  where  the  microfocus  x- 
ray  tube  was  operated  at  high  power,  80  W,  to  enlarge  the  diameter  of  the  x-ray  source  to 
approximately  75  pm.  With  this  source  diameter,  the  phase  contrast  component  of  the 
image  disappears  as  a  result  of  convolution  of  the  phase  contrast  fringes  over  the  source 
area  of  the  x-ray  beam.  The  exposure  time  for  the  image  was  75  s.  In  Fig.  5  (C)  a 
conventional  phase  contrast  image  of  the  two  beads  is  shown  using  the  same  total  x-ray 
fluence  as  was  used  in  Fig.  5  (B),  but  with  the  x-ray  tube  operating  at  10  W,  where  source 
diameter  was  25  pm.  The  fringes  at  the  perimeters  of  the  beads,  which  serve  to  define  the 
edges  of  the  beads,  are  the  result  of  interference  from  rapid  phase  variations  in  the  x-ray 
paths,  which,  for  a  sphere,  are  largest  at  its  perimeter 23 .  The  results  of  the  image 
subtraction  with  the  ultrasound  directed  onto  only  the  bead  closest  to  the  x-ray  tube  are 
shown  in  Fig  5  (D).  The  bead  that  was  not  irradiated  with  ultrasound  is  not  visible  in  the 
image,  while  the  phase  contrast  component  of  the  image  of  the  irradiated  bead  shows  up 
with  the  phase  contrast  component  highlighted.  The  subtracted  phase  contrast  image  in 
Fig.  5  (D)  shows  a  reversal  of  the  bright  and  dark  regions  at  the  perimeter  of  the  bead  in 
the  direction  of  the  ultrasound  at  the  two  sides  of  the  bead  (which  appear  at  the  top  and 
bottom  of  Fig.  5  (D)):  the  interior  of  the  arc  at  the  top  of  the  figure  appears  bright  while  it 


is  the  exterior  of  the  arc  at  the  bottom  of  the  figure  that  appears  bright.  The  reversal  of  the 
shading  is  a  consequence  of  subtraction  —  the  opposite  shading  could  be  produced  by 
reversal  of  the  order  of  subtraction  of  the  two  images. 


As  a  result  of  a  collaboration  with  the  Brown  Medical  School,  we  will  soon  be 
able  to  test  our  technique  on  biological  samples.  We  will  obtain  mouse  skin  and  liver 
tumors.  The  first  measurements  are  expected  to  corroborate  our  preliminary  proof  of 
concept  experiments. 


X-ray  beam 


Figure  5:  (A)  Side  view  of  the  sample  consisting  of  two  Teflon  beads  cast  in  agarose.  The  two  beads  are 
separated  by  3  cm.  (B)  Image  of  the  beads  with  no  phase  contrast.  (C)  Image  of  the  beads  with  phase 
contrast.  (D)  Subtracted  image  of  the  two  beads.  The  bead  irradiated  with  ultrasound  is  nearest  the  CCD 
camera  and  produces  a  slightly  larger  image  than  the  bead  located  outside  the  sound  beam.  The  x-ray  tube 
voltage  for  the  images  was  90  kV.  The  image  in  (B)  was  taken  with  a  75  s  exposure  time  with  the  x-ray 
tube  operating  at  80  W;  the  image  in  (C)  was  taken  with  eight  times  lower  tube  wattage  and  eight  times 
longer  exposure  to  give  the  same  total  x-ray  fluence. 
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X-ray  Phase  Contrast  Tomography 

Tomographic  back  projection  calculations  provide  the  reconstruction  of  a  2D 
object  from  a  serie  of  ID  projections 24 .  If  an  object  is  described  by  the  function  f(x,y), 
the  projection  p(0,r)  of  the  object  at  a  defined  angle  0  and  at  a  distance  r  can  be  written  as 
the  Radon  transform: 


p{6,r)=  Radon(f(x,y))=  J”  f(x,y)-S(x-cos6  +  ysin6-r)dxdy 

We  want  to  recreate  the  object  f(x,y)  knowing  the  projections  p(0,r),  namely  applying  an 
inverse  Radon  transform: 


f(x,  y)  =  Radon  1  ( Pr  oj(  x,  y)) 


2n  x-cosQ  +  ysinO -r 


If  all  the  projections  are  known,  one  could  theoretically  reconstruct  the  object.  However, 
the  inverse  Radon  equation  is  undefined  when: 
r  =  x-  cos  9  +  ysind 

Practically,  the  filtered  back  projection  algorithm  is  used  to  compute  the  inverse 
Radon  transform.  We  used  Matlab  7  to  simulate  and  reconstruct  defined  objects.  The 
filter  is  designed  directly  in  the  frequency  domain  and  then  multiplied  by  the  FFT  of  the 
projections.  The  projections  are  zero-padded  to  a  power  of  2  before  filtering  to  prevent 
spatial  domain  aliasing  and  to  speed  up  the  FFT. 

The  Radon  transform  can  be  used  to  describe  a  parallel  beam  geometry  as  shown 
in  Figure  6 A.  Our  experiment  approximates  a  point  source  with  a  fan  beam  geometry  as 
shown  in  Fig.  6B.  We  have  computed  our  simulation  transforming  our  fan  beam  data  into 
a  parallel  beam  geometry  using  Matlab  7  and  then  employing  the  Radon  and  inverse 
Radon  functions. 
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Parallel  Beam 


Fan  Beam 


Figure  7:  (A)  Numerically  simulated  phase  contrast  signal  of  a  hollow  disc  of  low  absorbing  material.  The 
expected  Fresnel  fringes  have  been  reduced  to  just  one  dark  and  one  light.  The  x-ray  source  is  on  the  left. 
The  axis  of  rotation  is  at  position  (500;500).  (B)  shows  the  fan  beam  projection  of  object  (A).  The  X  axis  is 
now  the  angle  value  and  the  Y  axis  is  the  projected  Y  axis  of  the  object  on  the  detection  plane.  (C)  shows 
the  back  projected  object. 

As  a  preliminary  test,  we  have  numerically  simulated  the  phase  contrast  signal  of 
a  well  defined  phase  object.  A  hollow  disc  of  low  absorbing  material  will  generate  a 
phase  contrast  signal  as  seen  in  Figure  7A.  This  represents  a  slice  of  a  3D  object  along 
the  Z  axis.  The  expected  Fresnel  fringes  have  been  reduced  to  just  one  dark  and  one  light 
ring.  The  x-ray  source  is  on  the  left  and  the  axis  of  rotation  is  at  position  (500;500).  The 
projection  shows  the  projected  intensities  of  the  object  on  the  detection  plane.  The 
scheme  is  as  follow,  from  a  2D  object,  one  obtains  ID  projections  at  defined  angles. 
These  ID  projections  are  put  together  in  Figure  7B.  This  matrix  is  then  used  to 
reconstruct  the  original  object  using  the  filtered  back  projection  algorithm,  Figure7C.  The 
reconstructed  object  is  clearly  in  agreement  with  the  original  projected  object. 
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Figure  8:  (A)  shows  a  projection  at  a  defined  high  and  angle  of  a  Teflon  tube.  (B)  shows  the  fan  beam 
projection  of  a  Teflon  tube.  The  X  axis  is  the  angle  value  and  the  Y  axis  is  the  distance  on  the  detector.  We 
measured  at  5  degree  increment.  (C)  shows  the  reconstructed  object  using  the  back  projection  algorithm. 

We  have  developed  an  experimental  setup  to  perform  3D  measurements  and  as  a 
first  test,  we  have  measured  a  standing  Teflon  tube.  This  sample  gives  a  comparable 
signal  to  the  simulated  hollow  disc  of  low  absorbing  material,  Fig.  7A.  Figure  8A  shows 
a  projection  at  a  defined  height  and  angle  of  the  Teflon  tube.  By  rotating  the  sample,  we 
obtained  the  projection  matrix,  Figure  8B.  Our  algorithm  reconstructs  2D  objects  from 
ID  projections.  In  the  case  of  a  3D  object,  the  typical  procedure  is  to  perform  the 
reconstruction  separately  at  different  heights  (z  axis)  and  put  the  2D  together  to  obtain  the 
final  reconstructed  3D  object.  Such  a  2D  reconstruction  is  shown  in  Figure  8C.  The  angle 
increment  of  5  degrees  is  responsible  for  some  discrimination  noise  on  the  finale  picture. 
We  can  now  combine  the  reconstructed  2D  objects  in  the  z  direction  and  obtain  a  series 
of  cross  sections  of  our  object. 

In  the  future,  our  tomographic  setup  will  be  used  in  combination  with  the 
acoustically  modulated  X-ray  phase  contrast  imaging  technique.  Some  interesting  3D 
pictures  of  ideal  objects  will  be  obtained;  later,  biological  samples  will  be  systematically 
investigated. 
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Key  Research  Accomplishments 

We  have  demonstrated  the  effects  of  an  acoustic  field  on  a  sphere  in  a  fluid. 
Furthermore,  we  have  performed  comparative  experimental  measurements  of  a  Teflon 
bead  in  water.  We  have  been  able  to  predict  successfully  the  radiation  force  on  the  Teflon 
bead  in  water. 

We  have  solidified  our  understanding  and  our  technique  by  performing  some 
additional  proof  of  concept  experiments.  That  is,  we  have  shown  that  an  object  is 
selectively  displaced  by  the  sound  field  and  not  by  secondary  collective  surroundings 
surface  or  volume  effects. 

We  have  demonstrated  the  ability  to  reconstruct  3D  object  using  filtered  back  projection 
algorithm. 

Reportable  Outcomes 

Best  Paper  Award  in  “Photons  plus  Ultrasound:  Imaging  and  Sensing”,  International 
Society  of  Optical  Engineering,  SPIE  given  by  Fairway  Medical  Technologies  and 
LaserSonix  Technologies,  San  Jose,  California,  2005,  to  be  published. 

A.  Beveridge,  C.  J.  Bailat,  T.  J.  Hamilton,  S.  Wang,  C.  Rose-Petruck,  G.  J.  Diebold  and 
V.  Gusev,  Imaging  with  Vibration  Potential  and  Phase-contrast  X-rays  (accepted  for 
publication) 

T.  J.  Hamilton,  C.  J.  Bailat,  C.  Rose-Petruck,  G.  J.  Diebold,  Acoustically  Modulated  X- 
ray  Phase  Contrast  Imaging,  Phys.  Med.  Biol.,  49  (2004)  4985-4996. 

C.  J.  Bailat,  T.  J.  Hamilton,  C.  Rose-Petruck,  G.  J.  Diebold,  Acoustic  radiation  pressure: 
a  phase  contrast  agent  for  x-ray  phase  contrast  imaging,  Appl.  Phys.  Lett.,  Vol.  85  No. 

19,  Nov.  2004,  4517-4519. 

Conclusions 

We  have  completed  all  the  necessary  construction  of  experimental  apparatus  and  proven 
the  concept  of  ultrasonically  modified  x-ray  phase  contrast  imaging.  We  have  shown  the 
method  has  the  ability  to  select  objects  within  a  body  for  imaging.  We  have  laid  the 
foundation  for  a  new  phase  contrast  tomography  with  ultrasonic  radiation  pressure. 
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